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L  INTRODUCTION 

In  the  early  1970s,  the  US  laser  fusion  program  was  centered  upon  a  simple  and 
elegant  pellet  concept:  a  small  bare  spherical  shell  of  deuterium-tritium  fuel  (DT) 
located  in  the  middle  of  a  large  empty  reactor  chamber,  directly  illuminated  by 
symmetrically  positioned  laser  beams  that  were  carefully  tuned  to  produce  high 
compression  and  central  ignition^  By  1975  this  concept  had  been  largely  abandoned. 
The  laser  beams  could  not  be  made  sufficiently  uniform  to  compress  the  shell 
symmetrically  inward;  the  target  design  required  very  high  laser  intensities  with  the 
certainty  of  triggering  numerous  laser-plasma  instabilities;  the  implosion  was  highly 
unstable  to  the  ablative  Rayleigh-Taylor  (RT)  hydrodynamic  instability;  and  the 
predicted  target  gain  was  too  low  to  overcome  the  low  efficiencies  of  realistic  lasers 
in  a  practical  reactor. 

Although  the  original  target  design  did  not  work,  the  general  concept  is  very 
appealing:  (a)  direct  illumination  of  a  spherical  target  maximizes  the  coupling 
efficiency  from  laser  energy  to  imploding  DT  fuel,  and  therefore  maximizes  the 
potential  target  gains;  (b)  tiie  spherical  divergence  of  the  plasma  corona  limits  the 
interaction  length  between  the  laser  light  and  tiie  plasma,  and  thus  diminishes  the 
risk  of  deleterious  laser-plasma  instabilities;  (c)  the  overall  spherical  symmetry 
minimizes  the  physical  complexity  of  the  target  and  thus  limits  the  number  of  ways 
that  it  can  fail;  (d)  the  pellet  mass  is  minimized,  which  then  minimizes  the  required 
laser  energy,  the  cost  of  the  pellet,  and  the  mass  of  the  debris. 

Through  the  years  there  has  been  a  continuing  low  level  of  effort  on  direct-drive 
laser  fusion,  and  now  the  above  physics  problems  appear  to  have  been  largely 
resolved.  The  laser  beams  can  be  optically  smoothed  to  produce  a  residual  rms  flux 
nonuniformity  that  is  well  below  one  percent.  The  risk  of  laser-plasma  instabilities 
has  been  minimized  by  lowering  the  peak  laser  intensity,  shifting  the  laser 
wavelength  into  the  UV,  and  using  optically  smoothed  laser  beams.  The  RT 
instability  can  be  ameliorated  by  preheating  the  ablator.  These  techniques  have  now 
been  pulled  together  to  g;enerate  three  specific  and  different  direct-drive  target 
designs  that  may  reach  ignition.  Some  may  have  high  enough  energy  gains  to  be 
useful  in  a  fusion  reactor.  These  designs  are  currently  being  evaluated  and 
optimized. 

This  paper  will  review  the  status  and  the  prospects  of  the  direct-drive  laser 
fusion  concept  as  of  November  1997.  We  will  summarize  the  experimental  and 
theoretical  advances  in  the  field,  describe  and  evaluate  flie  specific  target  designs, 
and  list  some  of  the  problems  that  remain  to  be  addressed. 

Although  construction  is  now  imderway  for  a  few-shots-per-day  megajoule  laser 
for  defense  applications,  there  has  been  very  little  development  of  lasers  that  are 
suitable  for  a  fusion  reactor.  These  fusion  lasers  would  need  a  much  higher 
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repetition  rate,  higher  durability,  higher  efficiency,  and  lower  cost.  Specific  laser 
concepts  have  been  proposed  that  may  meet  all  the  requirements,  using  either  a 
krypton-fluoride  laser  (KrF)  or  a  diode-pumped  solid-state  laser  (DPSSL). 

II.  LASER  BEAM  NONUNIFORMITIES 

A.  Laser  Beam  Uniformity  Requirements 

The  generic  direct-drive  pellet  is  shown  in  Fig.  1.  A  shell  of  frozen  DT  fuel  (a  few 
hundred  microns  thick)  is  surrounded  by  a  low-Z  ablator  that  is  burned  off  by  the 
laser  beams.  The  masses  of  the  DT  and  the  ablator  are  similar.  There  is  sometimes  a 
thin  outer  membrane  surroimding  the  ablator.  The  pulse  shape  has  two  temporal 
regions:  the  "foot"  and  the  "main"  or  "drive"  portions.  The  low-intensity  foot  of  the 
laser  pulse  generates  a  few  megabar  shock  wave  that  compresses  the  pellet  shell  and 
places  it  on  the  desired  isentrope.  The  laser  intensity  dien  rises  rapidly,  generating 
pressures  of  70-100  MBar  and  accelerating  the  shell  inward.  The  rising  pulse  shape 
is  carefully  timed  to  minimize  any  additional  increases  in  the  isentrope  of  tiie  DT 
fuel. 

Figure  2  shows  typical  density  and  temperature  profiles  at  final  stagnation,  witti 
die  low-density  central  spark  plug  surrounded  by  a  higher-density  main  fuel  region. 
The  hot  central  region  has  a  pR  ~  0.3  g/cm^,  approximately  equal  to  the  range  of  an 
alpha  particle  in  a  5  keV  spark  plug  plasma.  This  spark  plug  self-heats  because  of  the 
local  alpha  particle  deposition,  rising  to 


Fig.  1  Direct-drive  pellets  are 
geometrically  very  simple,  with 
approximately  equal  masses  in  the  DT 
fuel  and  the  ablator,  surrounding  a  DT 
gas  that  is  in  equilibrium  with  the  fuel. 
Sometimes  there  is  a  thin  outer  coating. 


Fig.  2  Typical  density  and  temperature 
profiles  for  a  direct  drive  pellet  at  final 
stagnation. 
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a  temperature  sufficient  to  initiate  a  propagating  bum  into  the  surrounding  colder 
fuel.  The  parameters  in  Fig.  2  were  taken  from  recent  calculations  of  a  direct-drive 
target  design  at  the  Naval  Research  Laboratory  (NRL).  The  spark  plug  radius  Rf  is 
~  4%  of  the  initial  pellet  radius.  If  one  postulates  that  any  angular  variation  SRf 
should  be  no  more  than  50%  of  Rf,  then  one  can  derive^  a  restriction  on  the  allowable 
low-mode  pressure  variations:  5P/P  <  0.5^Ry /Rq).  Since  P  ~  one  concludes  that 
the  very  low-mode  intensity  variations  around  the  pellet  should  be  less  than  2.5%. 

It  has  been  shown  analytically^  that  almost  any  even  number  (>  6)  of  carefully 
located  laser  beams  is  sufficient  to  produce  perfect  illumination  symmetry  on  a  hard 
sphere.  If  each  beam  generates  a  cos^  (6)  intensity  distribution  on  the  pellet,  then  the 
various  cos^  factors  can  add  to  unity.  When  one  includes  realistic  laser  profiles, 
refraction,  collisional  absorption,  inward  motion  of  the  pellet  surface,  laser  beam 
imbalances,  and  laser  beam  mispointing,  then  the  uniformity  can  no  longer  be 
perfect,  but  more  beams  will  improve  the  statistics.  Studies^  have  indicated  that  60 
or  more  laser  beams  should  meet  the  low  spherical  mode  {£  =  1-10)  symmetry 
requirements  for  an  implosion. 

Higher-mode  symmetry  criteria  are  more  severe,  because  the  high  density  shell  is 
being  accelerated  inward  by  a  lower  density  plasma  corona  and  the  interface  is  RT 
unstable.  For  the  simple  case  of  a  heavy  fluid  of  density  pz  supported  above  a  lower 
density  fluid  of  density  pi,  and  ignoring  finite  thickness  effects,  the  exponential 
growth  rate  of  an  interface  perhubation  is  given  hy  yo-  -yjAkg,  where  k  is  the  mode 
number  of  the  perturbation,  g  is  the  acceleration,  and  A  is  the  Atwood  number:  (pz- 
Pi)/ (P2+Pi)-  For  a  spherical  implosion,  the  spherical  harmonics  are  the  normal 
modes,  and  the  analogous  formula  for  large  £  can  be  approximated  by  substituting 
£  /R  for  k.  Thus  yo  =  ^A£g/R.  If  the  shell  accelerates  inwards  such  that  gt^  ~  Rq, 
then  the  number  of  e-folds  is  proportional  to  ^fA£.  The  higher  mode  numbers,  those 
produced  within  each  laser  beam,  have  the  larger  growth  rate.  The  modes  that  are 
most  dangerous  will  have  perturbation  wavelengths  that  are  comparable  to  the 
pellet  shell  thickness,  since  still  higher  mode  numbers  would  saturate  before 
breaking  through  the  pellet  shell.  To  reduce  the  mode  number  of  the  most 
dangerous  RT  modes,  target  designs  try  to  minimize  the  in-flight  aspect  ratio  R/ AR 
(average  shell  radius  divided  by  the  in-flight  shell  thickness). 

There  are  several  ways  of  reducing  the  growth  of  the  most  dangerous  modes, 
including  density  gradient  stabilization,  ablative  (convective)  stabilization,  nonlinear 
saturation,  mode-mode  coupling,  and  spherical  convergence;  these  will  be  described 
later  in  this  article.  Figure  3  plots  the  classical  number  of  RT  e-folds  J  Yo^t/  3Jod  the 
value  J  ydt  with  ablative  stabilization,  versus  mode  number  for  a  recent  target 
design.  The  most  unstable  modes  are  in  the  vicinity  of  100.  The  perturbation 
wavelength  [  A  =  27iR/£  =  2R{7t/£)]  is  ~n/100  of  the  size  of  an  individual  laser  beam. 
Since  this  wavelength  is  much  less  than  the  size  of  individual  laser  beams,  successful 
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fusion  requires  methods  of  reducing  the  transverse  structure  within  laser  beams  to  a 
low  level. 

B,  Optical  Smoothing  of  Laser  Beams 

1.  Idea  of  Beam  Smoothing  (Temporal  &  Spatial  Bandwidth) 


Large  laser  systems  have  many  optical  elements,  each  of  which  introduces  phase 
distortion.  At  the  focal  plane,  these  phase  variations  would  lead  to  interference 
effects  with  large  and  uncontrollable  illumination  nonuniformities.  The  solution  to 
the  nonuniformity  problem  was  based  on  the  observation  that  the  most  uniform 
illumination  is  not  coherent  laser  light,  but  rather  incoherent  light  from  conventional 
optical  sources  such  as  surdight.  Fusion  pellets  have  a  diameter  of  ~  5  mm,  while  the 
coherent  light  from  a  laser  can  ideally  be  focused  to  a  spot  size  of  oiUy  a  few 
microns.  We  can  therefore  trade  off  the  excess  coherency  and  focusability  of  a  laser 
to  achieve  better  uniformity. 

The  laser  beam  is  first  broken  up  into  multiple  independently-phased  beamlets 
whose  widths  d  are  small  compared  to  the  full  beam  width  D.  The  beam  divergence 
is  characterized  by  the  large  number  Ns  =  D/d  of  beamlets  across  the  beam.  For  a 
megajoule  class  fusion  driver.  Ns  >  100.  Each  beamlet  has  much  less  phase  distortion 
than  the  full  beam,  and  its  focal  profile  can  be  nearly  ideal. 

Because  of  interference  between  the  beamlets  in  each  beam,  the  instantaneous 


focal  spot  intensity  produces  tdghly 
nonuniform  random  speckle  whose 
envelope  is  the  ideal  distribution.  The 
spatial  frequencies  of  this  speckle, «  = 
(wjc/  Ky)  form  a  two-dimensional  array  of 
integer  multiples.  For  a  square  laser 
aperture,  the  ensemble-averaged  spatial 
power  spectrum  is  given  by 


1  10  100  1000 


mode  number 


where  InxI,  lnyl<Ns. 

Spatial  and  temporal  incoherence  is 
then  introduced  by  forcing  the  speckle 
patterns  to  change  randomly  on  a 
coherence  time  Tc  that  is  much  shorter 


Fig.  3  Typical  number  of 
exponentiations  of  the  ablative 
Rayleigh-Taylor  instability  during  an 
implosion.  The  upper  "classical"  curve 
assumes  a  growth  rate  proportional  to 

-\/to .  The  lower  curve  includes  ablative 
(convective)  stabilization.  The  maximum 
growth  is  near  mode  number  100. 
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than  the  effective  hydrodynamic  response  time.  In  the  ideal  case,  where  all  the 
beamlet  phases  change  independently,  all  of  the  speckle  modes  are  randomized  each 
coherence  time,  and  all  are  thus  smoothed  at  the  same  rate.  For  an  effective 
averaging  time  tav,  there  are  then  tav/ic  statistically-independent  patterns,  and  the 
above  equation  is  replaced  by 


a\n) 


(2) 


thus  the  speckle  nonuniformity  approaches  zero  at  sufficiently  long  averaging  times 
The  total  mean  square  nonuniformity  becomes  cr^  =  Tc/tav  The  effectiveness  of  the 
smoothing  process  thus  depends  strongly  on  the  laser's  bandwidth,  with  Av~l  /Tc. 
With  the  Nike  KrF  laser  the  bandwidth  has  reached  3  THz,  with  the  coherence  time 
as  short  as  0.35  ps. 


2.  Induced  Spatial  Incoherence 


The  most  effective  beam  smoothing  scheme  that  has  been  tested  so  far  is  the 
Induced  Spatial  Incoherence  (ISI)  technique^  A7A9.  a  spatially  and  temporally 
incoherent  light  source  uniformly  illuminates  the  front  end  (object),  which  is 
currently  chosen  to  be  a  hard  apq^re.  (This  aperture  could  also  be  apodized  to 
create  a  shaped  envelope  profile.)  The  object  profile  is  then  imaged  through  tiie  laser 
system  onto  the  target  to 
produce  the  focal  spot 
envelope.  To  avoid 
imprinting  gain 
nonuniformities  onto  the 
focal  spot,  each  amplifier 
stage  is  located  at  or  near 
the  object's  (and  target's) 

Fourier  transform  plane. 

ISI  smoothes  all  spatial 
frequencies  at  the  same 
rate,  as  described  in  Eq.  2. 


The  focal  spot  envelope 
will  remain  relatively 
insensitive  to  large  scale 
phase  aberration  imposed 
by  the  laser  system  if  the 
spot  size  is  many  times  the 
diffraction  limit,  and  also 


Fig.  4  The  measured  focal  spot  laser  intensity  is  very 
uniform  for  a  KrF  laser  beam  with  ISI  smoothing, 
with  1  THz  bandwidth,  and  with  a  4  ns  averaging 
time.  RMS  nonuniformities  in  each  beam  are  1%, 
excluding  tilt  and  curvature. 
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large  compared  to  the  point-spread  distribution  (PSD)  created  by  the  aberrations. 
The  Nike  KrF  laser  optical  design^®  allows  nearly  complete  compensation  of  all 
systematic  phase  aberrations,  such  as  astigmatism,  coma,  and  spherical  aberration. 

In  the  present  configuration,  the  PSD  width  due  to  random  phase  aberration  is  ~  15 
times  the  diffraction  limit.  For  the  75  XDL  focal  spots  currently  being  used,  this 
allows  a  flat  top  over  about  half  of  the  envelope,  as  shown  in  Fig  4.  The  flat  region 
could  be  extended  by  using  adaptive  optics  techniques,  such  as  those  used  in 
Beamletii  to  correct  for  the  phase  aberration  in  the  amplifiers. 

Using  ISI,  the  Nike  KrF  laser  has  produced  smooth,  flat-topped  focal 
distributions  with  residual  speckle  fluctuations  of  1%  rms  in  a  single  beam  with  1 
THz  bandwidth  and  Tc  =  0.66  ps  (for  a  Gaussian  spectrum,  s  0.66/Av).  With  37 
overlapped  Nike  beams,  the  nonuniformity  is  below  that  which  is  measurable,  but 
is  estimated  to  drop  by  as  much  as  1/V^  to  0.15%  rms  (ignoring  the  short 
wavelength  beam-beam  interference).  The  detailed  mode  spectrum  has  been 
measured  for  a  single  beani,  and  is  in  agreement  with  theory  12.  The  bandwidth  of 
Nike  for  some  recent  target  experiments  has  been  extended  to  3  THz. 

Equation  2  can  be  used  to  estimate  the  single-heam  contributions  to  the  low-mode 
nonuniformities  using  ISI  on  a  spherical  target.  For  modes  |M| «  N^,  Eq.  (2)  reduces 
to  the  azimuihally-symmetric  expression  (T^(M)  =  (t^  Th®  total  RMS 

nonuniformity  of  all  modes  up  to  some  value  |M|^  «  is  tihen 


If  the  target  diameter  equals  the  focal  beam  diameter,  dien  the  spherical  mode 
number  is  ^  =  7r|M],  and  the  corresponding  RMS  nonuniformity  of  all  -modes  up  to 
f  c  becomes 


,(|M|<|M|J  = 


Ns 


(3) 


^rms 


{e<Q= 


(4) 


As  an  example,  consider  the  illumination  geometry  of  the  National  Ignition  Facility 
(NTF),  where  a  3.2  mm  diameter  target  (diffraction  width  fX/ =  7 fjm  )  would 
require  Ns  =  457.  Without  any  temporal  smoothing,  the  nonuniformity  in  all  modes 
1  <  ^  <  28  would  then  be  approximately  28/(457V^j  =  3.5%.  If,  for  example,  one 
averages  over  100  ps  with  Tc  =  1  ps,  then  the  nommiformity  drops  to  0.35%.  Thus  the 
speckle  from  ISI  does  not  introduce  significant  low-mode  nonuniformities. 

Because  the  ISI  light  is  incoherent  and  chaotic,  its  speckle  varies  rapidly  in  space 
and  time  throughout  the  entire  laser  system.  This  is  not  a  problem  in  KrF  lasers. 
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where  the  gaseous  medium  and  low  saturation  flux  ensure  very  little  nonlinear 
phase  distortion  in  the  amplifiers.  However  it  would  be  undesirable  in  a  glass  laser, 
which  can  produce  large  nonlinear  phase  distortions,  and  which  also  requires 
harmonic  conversion. 

3.  Smoothing  hy  Spectral  Dispersion 

The  most  promising  beam  smoothing  scheme  for  solid-state  lasers  uses  instead  a 
technique  called  Smoothing  by  Spectral  Dispersion  (SSD)!^.  Her4  the  spatial 
divergence  is  created  near  the  final  focusing  lens  by  a  phase  plate  that  breaks  the 
beam  into  beamlets  with  adjustable  phase  shift.  The  temporal  modulation  is  created 
near  the  front  end  of  the  laser  by  a  phase  modulator  with  temporal  phase 
dependence  followed  by  a  diffraction  grating  that  spectrally  disperses  the  light 
over  an  angle  A6.  After  the  grating,  the  color  at  any  instant  cycles  across  the 
aperture,  and  the  interference  among  more  distant  beamlets  averages  out  in  time. 
Alternatively,  one  can  view  the  smoothing  process  as  the  result  of  moving  the 
speckle  pattern  back  and  forth  across  the  target  as  the  off-axis  angle  6(t)  changes  in 
time.  The  technique  has  been  extended  to  both  traitsverse  dimensions  (i.e.,  2-D  SSD) 
by  adding  a  second  modulator  (at  a  different  frequency)  and  an  orthogonal  grating 
to  disperse  the  light  in  the  second  transverse  direction. 

The  advantage  of  SSD  is  that  near  the  image  planes  of  the  gratings  the  variations 
are  just  in  the  phase,  thereby  avoiding  the  self-focusing  and  nonlinear  spectral 
broadening  that  could  limit  high  power  glass  laser  amplification.  Its  main 
disadvantage  so  far  has  been  that  it  cannot  effectively  smooth  the  lowest  spatial 
frequency  components  of  the  speckle  pattern,  which  arise  from  interference  among 
nearby  beamlets.  For  £  values  below  about  30,  the  mode  intensities  approach  the 
magnitudes  given  by  Eq.  (1)  rather  than  Eq.  (2). 

SSD  has  now  been  extended  to  two  transverse  dimensions  on  the  OMEGA  laser. 
In  addition  to  the  phase  plate,  there  is  also  a  polarization  rotatoria.  This  rotator  is  a 
birefringent  wedge  placed  before  the  random  phase  plate  to  vary  the  polarization 
across  the  beam,  reducing  the  nonuniformities  by  With  the  polarizer,  the 

single  beam  nonuniformity  is  expected  to  be  8.5%.  With  the  60  beam  of  OMEGA,  the 
statistical  overlap  should  reduce  the  nonuniformity  at  the  pellet  to  ~3.5%. 

Steps  are  being  taken  to  improve  the  degree  of  beam  smoothing  with  SSD. 
Techniques  have  been  proposed  to  increase  the  bandwidth  Av,  which  had  been 
limited  primarily  by  dispersion  in  the  sum-frequency  output  crystal.  One  can  extend 
the  phase  matching  bandwidth  of  this  crystal  by  replacing  it  by  a  pair  of  stagger- 
tuned  crystals^^. 

In  addition  to  a  broader  bandwidth,  the  spatial  phase  variation  across  the  beam 
(number  of  color  cycles)  must  be  increased  to  bring  the  smoothing  rate  of  the  longer 
wavelength  speckle  closer  to  the  ideal  rate  described  by  Eq.  2.  According  to  recent 


simulations^^,  with  increased  color  cycling  the  smoothing  should  approach  the  ideal 
rate  at  nearly  all  spatial  frequencies  down  to  the  limit  imposed  by  the  dispersion 
angle  A0. 

4.  Other  Smoothing  Techniques 

The  pure  ISI  technique  has  also  been  applied  to  the  Phebus  glass  laser  at  Limeil 
Laboratory^^.  The  laser  beam  uniformity  was  excellent,  and  agreed  with  ISI  theory, 
however  the  laser  had  to  operate  at  reduced  energy  to  avoid  nonlinear  optical 
distortion  in  the  amplifiers.  At  Osaka  University^^  an  ISI  beam  has  been  propagated 
through  the  Gekko  XII  laser  chain,  and  then  passed  through  a  phase  plate  to  increase 
the  beam  divergence,  producing  a  very  uniform  laser  beam  at  the  second  harmonic. 

III.  HYDRODYNAMIC  INSTABILITIES 

A.  Laser  Imprinting  and  the  Rayleigh-Taylor  Instability 

The  growth  of  the  mass  perturbations  for  each  spherical  harmonic  (ignoring  the 
m-variation)  can  be  written  as  Af  =  exp^J  X/dfj.  The  coefficient  has  four 

sources  from  the  pellet  and  laser:  (1)  perturbations  on  the  inner  surface  of  the  DT 
cryogenic  (ice)  layer;  (2)  perturbations  on  the  outer  surface  finish  of  the  capsule 
ablator  material;  (3)  long-wavelength  laser  perturbations  associated  with  beam-to- 
beam  power  imbalance,  mispointing,  and  misfocusing;  and  (4)  short-wavelengtii 
laser  perturbations  resulting  from  imperfections  in  individual-beam  spatial 
distributions.  There  is  also  nonlinear  mode  coupling. 

Success  with  direct-drive  requires  that  one  minimize  both  the  coefficient  A®  and 
the  growth  rate  y^.  Since  there  is  often  a  target  performance  penalty  associated  with 
designing  a  pellet  with  lower  y^,  the  goal  has  been  to  first  reduce  A°  to  as  low  a 
value  as  possible. 

The  laser  fusion  program  has  made  considerable  progress  in  the  fabrication  of 
smooth  pellets.  Recent  experiments^^  have  shown  that  it  should  be  possible  to  limit 
the  roughness  on  the  iimer  surface  of  the  DT  ice  to  ~  1.2  ±  0.2  pm.  Currently 
manufactured  pol5nner  shells^o  have  an  outside  surface  finish  whose  roughness  is 
dominated  by  the  lowest  modes.  For  mode  numbers  greater  than  ten,  the  rms 
roughness  is  in  the  range  of  90  A.  Some  pellet  designs  now  use  a  low-density  foam 
ablator.  The  surface  finish  of  foams  is  currently  ~  1  pm  rms,  but  some  further 
improvement  is  possible^!.  Although  full  success  in  target  fabrication  has  not  yet 
been  demonstrated  in  the  appropriate  geometries,  it  is  possible  that  target  roughness 
will  not  be  the  dominant  factor  in  the  coefficient  A®. 
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If  the  target  fabrication  can  be 
sufficiently  improved,  and  if  the 
lasers  can  be  aimed  with  sufficient 
accuracy,  then  the  dominant  part  of 
A®  will  arise  from  the  inevitable 
nonuniformities  within  each  of  the 
laser  beams.  Although  the  target 
or\ly  responds  to  a  smoother  time- 
average  of  the  laser  fluctuations,  this 
time-average  is  still  non  zero  in  any 
finite  time  interval-  The  laser- 
induced  mass  nonuniformity^  is 
called  "laser  imprinting,"  and  is 
especially  large  in  the  earliest 
portion  of  the  laser  pulse,  when  there 
is  only  a  thin  plaisma  corona 
separating  the  laser  absorption 
surface  from  the  ablation  surface. 
Later  in  time,  thermal  conduction 
and  acoustic  waves  in  the  thicker 
corona  reduce  the  laser-imposed 
nonuniformities . 

Figures  5a  and  5b  show  the  laser 
imprinting  versus  time  from  a  two- 
dimensional  (2D)  computer 
simulation  of  a  pellet  implosion  that 
includes  a  model  of  the  time- 
dependent  ISI  laser  fluctuations.  The 
higher  modes  saturate  or  oscillate 
early  in  time,  while  the  lower  modes 
continue  to  grow  during  the  entire 
foot  of  the  pulse,  which  is  ~  20  ns  in 
this  case. 

During  the  high  intensity  portion 
of  the  laser  pulse  the  shell  accelerates 
inward  and  the  most  dangerous 
perturbation  modes  grow 
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Fig.  5  T)pical  growth  of  mass 
perturbations  in  various  spherical 
harmonics  during  a  pellet  implosion. 
During  the  first  20  ns  this  laser  pulse  has  a 
low  laser  intensity  foot,  followed  by  a  high 
intensity  main  pulse.  During  the  foot  the 
laser  "imprints"  perturbations.  The  lower 
mode  numbers,  32, 64  and  96  grow  during 
the  entire  foot,  and  then  rise  due  to  the 
Rayleigh-Taylor  instability  in  the  main 
pulse.  The  higher  modes,  >  128,  are 
saturated  by  mass  ablation  during  the  foot, 
and  then  grow  during  the  main  pulse. 


exponentially  due  to  the  RT  instability.  The  largest  amplitude  perturbations  have 


wavelengths  comparable  to  the  in-flight  shell  thickness,  with  mode  numbers  in  the 
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range  of  50-100.  Saturation  of  these  modes  occurs  either  because  the  spherical 
convergence  closes  up  the  bubbles,  or  because  the  target  is  destroyed. 

The  early  time  imprinting  on  the  pellet  surface  is  equivalent  to  a  driven 
Richtmyer-Meshkov  instability,  with  saturation  due  to  mass  ablation^s.  Several 
methods  are  under  consideration  to  reduce  this  laser  imprinting.  (1)  Increasing  the 
laser's  bandwidth  should  reduce  the  imprint  in  any  given  time  interval  by  ~ 

1/  VA<d.  (2)  Lowering  the  density  of  the  ablator  reduces  the  imprint^^  by  (this  is 

somewhat  balanced  by  an  increase  in  the  foot  duration  for  the  thicker  target).  (3) 
Using  a  longer  wavelength  laser^^  during  the  foot  of  the  pulse  increases  the 
separation  between  the  laser  absorption  region  and  the  dense  ablator.  (4)  Adding  a 
"foam  buffer"  and  a  high-Z  radiator  to  the  outside  of  the  ablator^  creates  a  smooth 
supersonic  x-radiation  pulse  that  produces  an  early-time  plasma  corona  buffer.  (5) 
Adding  a  high-Z  coating  directly  to  the  surface  of  the  ablator  adds  extra  free 
electrons,  which  also  increases  the  early-time  separation  between  die  absorption  and 
the  ablation  regions.  The  application  of  some  of  these  techniques  to  specific  target 
designs  is  described  in  Section  DI. 

In  addition  to  a  minimization  of  laser  imprinting,  success  witii  laser  fusion 
requires  that  the  RT  growth  be  reduced  to  an  acceptable  level.  There  are  two  basic 
mechanisms  to  reduce  the  growth  rate:  gentle  density  gradients  and  convective  flow. 

If  the  minimum  density  gradient  scale  length  =  min|p/(dp/dx)[  and  the 
perturbation  wave  number  k  satisfy  kL^^  » 1,  then  the  growth  rate  is  reduced26/27,28 
to  ^g/L^ .  Radiation  can  enhance  L^. 

In  laser  fusion  the  ablation  process  convects  material  through  the  unstable 
region.  If  matter  can  flow  the  distance  of  one  perturbation  wavelength  within  one  e- 
folding  time  then  the  perturbations  are  swept  out  of  the  unstable  region  and  are 
convectively  stabilized.  Thus  RT  stabilization  should  occur  i(  v^fYg>l/k,  where  Va 
is  mass  ablation  velocity.  One  can  therefore  expect  that  there  would  be  an 
approximate  RT  equation  for  the  growth  rate  of  the  form 

r  =  oc^^r^|j^-pkv, 

The  ablative  stabilization  concept  and  the  above  equation  wais  first  derived 
analytically 29  from  a  crude  model  that  predicted,  pessimistically,  that  =  1.  In  recent 
years  the  commimity  has  used  a  better  numerical  fit30,31,32,33  34  vvith  =  0.9  and 
p  =  3.  However  this  newer  value  of  is  only  applicable  when  radiation  can  be 
neglected^^  xhe  theory  of  ablative  stabilization  of  RT  has  now  been  extensively 
studied  by  several  groups  35,36373839,40,41,42,43,44,45  various  detailed  formulas 
have  been  derived.  The  newer  expressions  can  still  be  reasonably  approximated  by 
the  above  equation  for  y,  with  the  coefficients  and  now  dependent  upon  the 
target  composition.  Reference  46  provides  a  review  of  the  various  models  of  ablative 
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stabilization,  along  with  useful  estimates  for  the  coefficients  and  in  different 
regimes. 

The  RT  growth  can  be  reduced  by  increasing  the  velocity  Va  at  the  ablation 
surface,  in  the  second  term  of  the  equation.  This  velocity  is  related  to  the  density  pa 
at  the  ablation  front  through  the  equation  for  the  mass  ablation  rate  dm/dt  =  pa  Va- 
To  enhance  Va  one  therefore  reduces  the  density  pa,  since  dm/dt  is  relatively 
insensitive  to  the  details  of  a  target  design.  For  a  given  ablation  pressure,  the  density 
is  reduced  when  the  target  is  at  a  higher  temperature.  Because  both  the  density  and 
temperature  of  a  target  vary  continuously  during  the  laser  pulse,  a  better  parameter 
to  use  is  the  isentrope,  proportional  to  P/p’'.  Because  the  Fermi-degenerate  pressure 
varies  as  p^,  it  is  convenient  to  use  the  iseritrope  parameter  a  =  P/Ppp,  where  P^p  is 
the  Fermi-degenerate  electron  pressure  of  a  fully  ionized  system.  (Note  that  this 
isentrope  a  is  not  the  same  as  the  cXpp  in  the  RT  equation.) 

Designing  targets  to  suppress  the  RT  instability  is  then  equivalent  to  designing 
targets  with  a  larger  a  at  the  ablation  front.  However  raising  a  can  also  reduce  the 
compressibility  and  the  fractional  bum  up  of  the  DT  fuel.  For  the  best  fuel  • 
compression  one  wants  to  keep  a  near  one.  Some  targets  are  proposed  here  that 
have  a  non  uniform  spatial  profile  for  a,  maximized  in  the  ablator  region  md 
minimized  in  the  DT  fuel. 

There  are  two  basic  methods  for  adjusting  the  isentrope  of  an  ablator:  radiation 
preheating  and  shock  preheating.  Applications  of  these  techniques  to  specific  target 
designs  are  described  in  Section  IE. 

In  the  first  stage  of  target  design  we  use  a  ID  spherical  radiation/hydrodynamics 
code,  along  with  post-processed  analytic  dispersion  relations  to  estimate  the  RT 
growth.  These  RT  relations  usually  include  initial  imprinting,  density  gradients, 
ablative  stabilization,  and  saturation^7'48,49,50,5i.  However  these  RT  relations  are  only 
approximate,  and  two-dimensional  (2D)  or  three-dimensional  (3D)  implosion  codes 
are  needed  for  more  accurate  analysis.  The  2D  computer  models  we  use  have  a  wide 
variety  of  physics  packages,  including  the  spatial  and  temporal  incoherence  of  the 
laser  light  with  specified  lens  sources,  opacities  and  emissivities  of  x-radiation  from 
atoms  that  are  not  in  thermodynamic  equilibrium  (non-LTE),  detailed  equations  of 
state  for  both  low  md  high  density  materials,  nonlinear  evolution  of  RT  into  bubbles 
and  spikes,  and  spherical  convergence. 

B.  Experimental  evaluation  of  laser  imprinting  and  Rayleigh-Taylor 

Because  of  the  inherent  limitations  in  both  the  numerics  and  the  physics  packages  of 
any  computer  modeling  there  is  of  course  a  need  for  experiments,  especially  to  test 
the  modeling  of  laser-imprinting  and  RT  growth.  Direct-drive  experiments  to 
measure  the  RT  growth  have  been  performed  for  a  dozen  years  at  various 
labs52^53^54  55,  However  the  diagnostics  in  these  experiments  could  not  reach  down 
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into  the  linear  perturbation  regime,  and  the  measurements  were  strongly  influenced 
by  the  relatively  large  nonimiformities  in  the  driving  lasers  and  the  lack  of  pulse 
shaping.  Only  very  recently  have  experiments  been  able  to  combine:  (a)  improved 
two-dimensional  laser  beam  smoothing  using  ISI  or  2-D  SSD;  (b)  few-nanosecond 
laser  pulse  lengths  with  pulse  shaping  to  simulate  the  hydrodynamics  of  a  fusion 
target;  and  (c)  improved  spatial  resolution  that  can  follow  the  shorter  wavelength 
perturbations;  however  we  have  still  not  reached  down  sufficiently  into  the  linear 
perturbation  regime.  These  newer  experiments,  summarized  below,  were  performed 
over  the  past  two  years  with  the  OMEGA  glass  laser  at  the  University  of  Rochester 
Laboratory  for  Laser  Energetics  (URLLE),  and  the  Nike  KrF  laser  at  the  Naval 
Research  Laboratory  (NRL). 

Since  the  most  dangerous  RT  wavelengths  are  comparable  to  the  pellet  shell 
thickness,  and  small  compared  to  the  shell  radius,  one  can  reasonably  approximate 
the  early-time  growth  of  the  perturbations  with  flat  foil  targets.  To  minimize  edge 
effects,  these  flat  foils  are  only  accelerated  a  distance  up  to  approximately  one-half 
the  laser's  spot  size.  Flat  foil  targets  also  have  the  advantage  of  easier  diagnostic 
access.  By  creating  a  separate  backlighter  that  propagates  x-rays  through  the  foil, 
one  can  directly  measure  the  time-dependent  transverse  mass  variations  in  the 
target.  By  using  a  separate  side-on  x-ray  source,  one  can  also  directly  measiu-e  the 
time-dependent  acceleration  of  the  target. 

1.  Evaluation  with  modulated  targets 

By  first  imposing  a  single-mode  sinusoidal  perturbation  and  comparing  its 
growth  with  code  predictions,  one  can  evaluate  whether  these  codes  can  correctly 
model  the  complex  hydrodynamic  behavior,  and  whether  the  laser  conditions  are  as 
expected. 

The  recent  OMEGA  data  are  shown  in  Fig.  6.  These  experiments  used  a  3  ns 
square  pulse  with  a  peak  intensity  of  2x10^^  W/cm^  on  plastic  (CH)  targets  that 
were  20  pm  thick,  with  a  1000  A  aluminum  overcoat  on  the  laser  illiunination  side  to 
keep  the  material  opaque  to  the  beginning  of  the  laser  pulse.  Data  are  shown  for 
sinusoidal  wavelength  perturbations  of:  a)  60  pm,  b)  31  pm,  and  c)  20  pm.  The  initial 
amplitudes  of  these  sinusoidal  modulations  were:  a)  0.5  pm,  b)  0.45  pm,  and  c)  0.5 
pm.  The  experimental  data  (plotted  as  diamonds)  and  the  ORCHID^^  predictions 
(plotted  as  a  line)  are  in  reasonably  good  agreement  for  all  wavelengths.  The 
rollover  of  the  measured  perturbation  optical  depth  at  late  times  has  been  shown  by 
the  simulations  to  be  due  to  the  inability  of  the  diagnostic  to  resolve  the  finer  spikes 
of  dense  material  at  the  ablation  surface.  The  increased  scatter  in  the  data  for  a  20  pm 
wavelength  perturbation  is  due  to  the  resolution  limits  of  the  x-ray  framing  camera. 
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Similar  experiments  have  been 
performed  with  the  Nike  laser  using  a 
somewhat  lower  laser  intensity  of 
W/cm^,  but  with  a  multi¬ 
nanosecond  foot  to  keep  the  target  on  a 
lower  isentrope  The  2D  FAST 
computer  simulations^^  were  also  in 
reasonably  good  agreement  with  the 
measured  longer  wavelength  initial 
perturbations,  but  deviated  from  the 
shorter  wavelength  perturbations  when 
the  long  RT  spikes  are  expected  to  break 
into  a  three-dimensional  structure  that 
can  not  be  modeled  with  a  2D  code. 

The  relatively  good  agreements  of 
OMEGA  and  Nike  experiments  with  the 
2D  modeling  suggests  that  these  codes 
can  simulate  the  basic  thermal  transport 
and  hydrodynamic  response  of  direct- 
drive,  finite-amplitude,  longer- 
wavelength  initial  modulations.  There 
are  however  three  potential  weaknesses 
in  the  above  experiments.  The  peak  laser 
intensity  was  only  2xl0l‘^  W/cm^, 
whereas  the  current  fusion  targets  have 
a  higher  peak  intensity,  ~  5-10  xlO^^ 
W/cm^.  At  lower  intensities  the  heat 
transport  can  be  characterized  by  a 
Spitzer-Harm  type  thermal 
conductivity^^.  At  higher  laser 
intensities  and  higher  electron 
temperatures  the  tail  of  the  Maxwellian 
hot  electrons  can  have  a  long  enough 
mean  free  path  to  directly  preheat  the 
target,  thereby  reducing  the  RT  growth 
rate.  Fokker-Planck  computer  modeling 
and  RT  experiments  with  short 
wavelength  lasers  have  not  yet  fully 
resolved  this  uncertainty.  Future 
experiments  at  URLLE  can  study  the  RT 
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Fig.  6  Evolution  of  the  optical  depth  as 
a  function  of  time  from  an  initial 
perturbation  on  the  ablation  surface  for 
planar  Rayleigh-Taylor  experiments 
using  the  OMEGA  laser  with  a  3  ns 
square  laser  pulse.  Data  are  shown  for 

(a)  60  pm,  (b)  31  pm,  and  (c)  20  pm 
wavelength  perturbations.  The 
diamonds  are  the  experimental  data  and 
the  solid  lines  are  the  ORCHID 
numerical  simulations.  The  square  point 
at  time  =  0  is  the  initial  amplitude 
multiplied  by  the  x-ray  mass  absorption. 
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behavior  at  higher  intensities  by  imploding  spherical  targets. 

The  second  potential  weakness  in  the  RT  experiments  has  been  the  choice  of  solid 
CH  as  the  target  material.  This  was  convenient  for  initial  experiments,  and  provided 
a  useful  comparison  of  our  experimental  and  computational  capabilities.  Future 
experiments  will  shift  to  the  use  of  the  ablator  materials  to  be  used  in  fusion  targets: 
cryogenic  DD  and  DT,  DD  wicked  into  CH  foams,  empty  CH  foams,  high-Z 
coatings,  etc. 

The  third  potential  weakness  is  the 
use  of  2D  codes  to  simulate  a  3D 


phenomenon.  The  experimental 
saturation  level  of  the  RT  mass 
modulations  is  less  than  the  predictions 
from  our  2D  codes.  This  should  not 
provide  a  basis  for  optimism.  Most 
likely,  the  long  thin  RT  spikes  break  up 
and  flow  into  the  bubble  regions, 
reducing  the  measured  mass 
nonuniformities.  Simulations  with 
NRL's  3D  code^'  predict  that  die 
saturation  level  of  the  ablative  RT 
instability  can  be  as  much  as  40%  higher 
than  predicted  with  2D  codes.  The 
predictions  from  this  3D  code  should  be 
taken  seriously  because  they  have  been 
successfully  compared^  with  laser 
imprinting  experiments  on  the 
VULCAN  laser.  Until  3D  codes  are  more 
widely  used  to  design  target  implosions, 
a  rule  of  thumb  should  probably  be  that 
2D  target  designs  include  a  margin  of 
safety  to  accoimt  for  RT  effects. 

2.  Evaluation  of  laser  imprinting  with 
smooth  targets 

The  next  stage  of  the  experimental 
program  has  been  to  accelerate  smooth 
targets  without  initial  mass 
modulations,  to  determine  if  the  3D 
laser  imprinting  of  incoherent  light  from 
real  lasers,  and  the  subsequent  RT 


time  (ns) 

Fig.  7  Measured  Rayleigh-Taylor 
perturbations  from  acceleration  of  a 
smooth  planar  target.  The  perturbations 
arose  from  imprinting  by  the  Nike  laser. 
The  squares  are  the  experimental  data 
and  the  two  solid  lines  are  the  FAST 
numerical  simulations.  The  dashed  line 
is  for  a  target  with  100  A  surface  finish, 
and  a  perfectly  uniform  laser.  The  solid 
line  is  for  a  perfectly  uniform  surface 
finish,  witii  37 ISI  laser  beams.  The  laser 
pulse  was  shaped  to  keep  the  target  on  a 
lower  isentrope.  The  use  of  a 
monochromatic  backlighter  allows  the 
data  to  be  converted  from  optical  depth 
to  mg/cm^. 
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growth,  is  adequately  modeled  by  the  2D  codes. 

With  ultra-uniform  laser  beams,  the  laser  imprinting  in  the  foot  of  the  pulse 
produces  a  mass  nonuniformity  that  is  too  small  to  directly  measure,  typically 
equivalent  to  100-1000  A  of  surface  modulation.  The  imprinting  is  only  observable 
after  being  amplified  by  the  RT  instability  during  the  main  portion  of  the  laser  pulse. 
In  these  experiments  the  RT  instability  is  therefore  used  primarily  as  a  diagnostic  of 
the  imprinting. 

Figure  7  shows  some  measurements  of  the  RT-amplified  imprinting  using  37 
overlapped  Nike  laser  beams.  These  CH  targets  were  40  |im  thick,  with  a  very  small 
initial  surface  roughness  (20-30  A  rms).  A  3-4  ns  foot  on  the  Nike  pulse  with 
amplitude  of  2%  of  the  main  pulse  simulated  the  long  low  intensity  foot  of  high  gain 
target  designs.  This  foot  was  followed  by  a  4  ns  main  pulse  that  accelerated  the 
target  at  0.8x10^^  W/cm^.  The  foot  and  main  pulses  were  chosen  such  that  the 
shocks  arrive  almost  simultaneously  at  the  rear  side  of  the  target,  thereby  keeping 
the  target  on  a  low  isentrope  (a  is  predicted  to  slowly  increase  from  2  to  4  during  the 
main  pulse)  .The  experiment  employed  a  spherically-curved  crystal  imaging  optic^ 
with  the  x-ray  backlighter.  This  technique  produces  a  nearly  monochromatic  image 
of  the  backlit  target  (at  1.86  keV)  so  that  the  optical  depth  modulations  could  be 
converted  into  absolute  areal  mass  nonuniformities. 

Figure  7  shows  the  areal  mass  nonuniformities  from  three  laser  shots  compared 
with  two  FAST  2D  simulations:  (a)  ISI  laser  imprinting  with  a  perfect  target;  and  (b) 
perfect  laser  beams  with  a  target  that  had  a  100  A  rms  initial  surface  roughness  (laser 
illumination  side).  There  is  relatively  good  agreement  between  experiment  and  the 
ISI  simulation  modeling  for  the  time  of  appearance  of  measurable  mass 
nonuniformities  (more  evident  if  Fig.  7  were  replotted  on  a  log  scale).  The  saturation 
amplitude  is  however  2-3  times  smaller  than  predicted.  This  difference  may  be  due 
to  3D  effects  at  the  shorter  transverse  wavelengths  that  can  not  be  modeled  by  a  2D 
simulation.  The  results  also  indicate  that  the  laser  imprinting  can  be  very  small  with 
state-of-the-art  beam  smoothing,  equivalent  to  about  100  A  rms  initial  surface 
roughness  for  these  Nike  parameters. 

III.  IGNITION  AND  HIGH-GAIN  TARGET  CONCEPTS 

Figure  8  illustrates  three  different  pellet  concepts  for  direct-drive  laser  fusion. 

The  first  design  using  a  DT  ablator  was  proposed  by  C.  Verdon  of  URLLE  at  the 
November  1993  APS  Plasma  Physics  meeting^^  and  is  published  here  for  the  first 
time.  The  intensity  of  the  first  shock  has  been  raised  to  put  the  DT  on  a  higher 
isentrope.  Typically,  the  DT  has  an  isentrope  a  ~  2-3.  This  target  is  calculated  to 
have  acceptably  low  laser  imprinting  and  RT  growth,  with  an  energy  gain  of  20 
using  a  1.5  MJ  laser  such  as  the  NBF.  This  target  has  two  major  advantages.  First, 
because  the  ablator  and  the  fuel  consist  of  the  same  material,  there  is  no  risk  of  a 
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classical  RT  instability  at  the  ablator/ fuel  interface.  Second,  making  the  ablator  out 
of  DT  maximizes  the  rocket  efficiency  of  the  implosion,  and  thus  maximizes  the 
potential  energy  gains.  There  is  however  one  basic  disadvantage.  Because  the  shock 
preheating  of  the  ablator  also  raises  the  isentrope  of  the  fuel,  there  is  reduced  fuel 
compression  and  thus  reduced  energy  gains,  compared  to  keeping  the  DT  fuel  on  a 
lower  isentrope. 

In  the  second  target  the  ablator  has  been  replaced  with  a  low  density  CH  foam. 
The  foam  is  again  preheated  with  an  initial  shock,  in  a  fashion  similar  to  the 
previous  target.  However  the  ablator  can  now  be  put  on  a  higher  isentrope  than  the 
DT  fuel.  With  this  spatial  profile  modification  of  a,  one  can  hopefully  enhance  RT 
control  while  at  the  same  time  minimizing  the  preheat  of  the  fuel.  This  target  has 
been  proposed  and  is  being  evaluated  by  A.  Schmitt  and  J.  Gardner  of  NRL.  It  is 
described  in  greater  detail  in  Ref.  51.  Implosion  studies  with  FAST2D  indicate  that 
the  combined  laser  imprinting  and  RT  growth  are  acceptable,  with  tcuget  gains  of 
approximately  80  for  a  3.5  MJ  laser.  The  potential  disadvantage  of  this  design  has 
been  its  sensitivity  to  the  outer  surface  finish  of  the  low-density  foam  ablator. 

The  ablator  in  the  third  target  is  DT  wicked  into  a  low-density  CH  foam,  along 
with  a  coating  of  a  high-Z  material  such  as  tungsten.  Instead  of  shock  preheating,  it 
uses  radiation  preheating.  During  the  low  intensity  of  the  foot  pulse,  tungsten  emits 
strongly  at  200-250  eV,  below  the  K-edge  of  carbon.  This  radiation  penetrates  deeply 
into  the  low-opacity  ablator.  This  target  has  die  advantage  that  the  fuel  can  be  kept 
on  a  very  low  isentrope  while  the  ablator  is  on  a  high  isentrope.  The  tungsten  blows 
off  during  the  foot  of  the  laser  pulse,  and  the  main  laser  pulse  directly  illuminates 
the  mostly-DT  ablator,  thus  maintaining  the  advantages  of  the  first  target.  Because 
the  mass  of  the  ablator  is  mostly  DT,  it  is  less  sensitive  to  the  surface  finish  of  die 
foam.  This  third  design  is  proposed  by  S.  Bodner  and  J.  Gardner  of  NRL  and  is  now 

foam 
ice 

DT 
vapor 


DT  ablator,  Foam  ablator, 

shock  preheat  shock  preheat 

Fig.  8  Three  possible  direct  drive  target  designs  that  may  reach  ignition,  (a)  An  all 
DT  design  with  shock  preheating;  (b)  A  low  density  CH  foam  ablator  with  shock 
preheating;  (c)  An  ablator  of  CH  foam  plus  frozen  DT,  preheated  by  radiation  from  a 
high-Z  coating. 


high  Z  coating 


DT  +  Foam  ablator, 
radiation  preheat 
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being  more  thoroughly  evaluated.  The  potential  disadvantages  are  not  yet  known. 


A.  Target  design  with  shock-preheated  DT  ablator 

1.  Simplified  Model  to  Identify  Acceptable  Design  Space 

A  very  simplified  constant-acceleration  model  can  be  used  to  quickly  estimate, 
for  a  given  RT  instability  model,  the  amount  of  perturbation  growth  a  direct-drive 
capsule  design  will  experience  during  the  acceleration  phase  and  how  this  growth 
compares  to  the  over-dense  shell  thickness  during  the  implosion. 

As  discussed  above,  the  parameters  urt  and  fi  in  the  dispersion  relation  for  the 
RT  growth  rate  are  dependent  on  the  material  composition  of  the  ablator.  For  a 
nominal  density  CH  ablator  (p  =  1.07  g/cm^),  detailed  simulations  have  shown  that 
typically  =  0.98  and  fiva  =  2.6x105  cm/s.  For  a  DT  ablator  (p  =  0.25  g/cm^), 

Ogj.  =  0.95  and  fiva  =  8.1x105  cm/ s.  (These  values  are  almost  independent  of  laser 
intensity,  as  shown  in  Ref.  46.)  The  higher  value  for  fioa  is  due  to  a  higher  ablation 
velocity  and  a  steeper  scale  length. 

For  direct-drive  implosions  of  solid 
density  material,  DT  is  the  optimum 
ablator  for  suppression  of  the  RT 
ablation  surface  instability.  In 
addition,  simple  scaling  arguments 
and  detailed  numerical  simulations 
have  foimd  that  for  directly  driven 
implosions  the  ablation  velocity 
scales  as  ««  where  a  is  the 
isentrope.  Using  this  information, 
along  with  the  saturation  model  of 
Ref.  44,  a  simple  model  can  be 
developed  for  the  tolerable  level  of 
initial  nonuniformity. 

This  model  assumes  a  constant 
acceleration  from  the  initial  radius 
Rq  to  Rq/ 4,  (at  which  point  the 
imploding  shell  begins  to  thicken 
due  to  spherical  convergence),  and  a 
spectrum  of  the  initial 
nonuniformity,  assumed  for  this 
example  to  be  of  the  form  C,{t)  = 
where  ^(£)  is  the  perturbation 
amplitude  for  the  Legendre  mode  £. 
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Fig.  9  Solid  curves:  inverse  mix-region 
width  of  a  pellet  [expressed  in  terms  of  an 

in-flight  ratio  (^/4)/Arnyx  ]  as  a  function  of 

degeneracy  parameter  a  for  a  number  of 
initial  perturbation  amplitudes.  Dashed 
curves:  in-flight  aspect  ratio  (^/4)/A^lieii 
as  a  function  of  a  for  capsule  designs  with 
two  different  implosion  velocities  Vjmp* 
a  given  a,  the  initial  perturbation  amplitude 
should  be  small  enough  that  the  operating 
point  on  the  solid  curve  lies  above  the 
dashed  curve  (i.e.,  <  Arsheii). 


The  model  calculates  die  width  of  the  Rayleigh-Taylor  mix  region  as  a 

function  of  time. 

Results  are  shown  in  Fig.  9.  The  solid  curves  give  (Ro/4)/dr^j.fas  a  function  of 
the  isentrope  a  for  various  initial  perturbation  amplitudes  Q  from  100  A  to  5000  A. 
For  a  given  decreases  as  the  departure  from  Fermi  degeneracy  (a)  increases, 

since  the  increased  ablation  velocity  leads  to  reduced  RT  growth.  The  dashed  curves 
give  (Ro/4) / Ar shell  for  two  capsule  implosion  velocities  close  to  3  x  10^  cm/ s. 

For  a  given  the  intersections  of  the  solid  curves  with  the  dashed  curve  give  the 
maximum  permissible  as  a  function  of  cl 

As  an  illustration,  a  capsule  design  with  an  isentrope  a  =  1  in  the  fuel  and 
ablator,  and  an  incident  energy  of  1.5  MJ,  produces  an  implosion  velocity  between 
2.5  and  3.5  x  10^  cm/s.  Figure  9  shows  that  the  maximum  perturbation  amplitude 
must  be  less  than  ~  100  A  if  the  mix  thickness  {Ar^i;^  is  not  to  exceed  the  over  dense 
shell  thickness  (Argjjgn)  during  the  implosion.  This  amplitude  represents  tibe 
combination  of  capsule  surface  finish  and  initial  laser  imprinting  added  in 
quadrature.  As  the  isentrope  a  is  increased,  the  amoimt  of  perturbation  that  the 
capsule  can  withstand  increases.  For  the  incident  energy  range  of  approximately  1.5 
MJ,  ignition  and  moderate  gain  require  an  a  no  greater  than  3  to  4.  From  Fig.  9  the 
maximum  initial  pierturbation  amplitude  must  then  be  less  than  5000  A. 

Detailed  two-dimensional  stability  calculations  have  shown  that  the  results  from 
this  simple  model  are  reasonably  accurate  for  the  acceleration  phase  of  the 
implosion.  When  the  earlier  portions  of  the  capsule  implosion  are  accurately  treated, 
the  initial  perturbations  must  be  reduced  by  a  factor  of  about  2  to  3  from  those 
shown  in  Fig.  9,  because  there  is  a  small  additional  amount  of  continual  growtfi  of 
longer-wavelength  perturbations  during  the  foot  of  the  pulse.  When  the  Rayleigh- 
Taylor  growth  of  initial  inside-surface  perturbations  during  the  deceleration  phase  is 
also  included,  the  acceptable  initial  outside-surface  perturbation  must  be  reduced  by 
another  factor  of  approximately  2.  Thus,  the  initial  perturbation  amplitude  must  be 
<1000  A  when  a =3. 

While  obviously  not  a  replacement  for  actual  multidimensional  simulations,  this 
model  has  been  found  to  be  an  accurate  indicator  of  whether  a  single-shell  direct- 
drive  capsule  is  likely  to  experience  excessive  shell  distortions  during  the  implosion 
phase  and/or  what  the  level  of  "initial"  perturbations  must  be  limited  to  in  order  to 
avoid  shell  disruption  during  the  implosion. 

2.  Ignition  and  High-Gain  Single-Shell  Designs 

In  this  section  we  will  consider  a  series  of  "all-DT,"  direct-drive,  single-shell 
cryogenic  capsules  designed  for  the  incident  energy  range  expected  for  the  NIF;  1.5 
to  1.8  MJ  at  351  nm.  These  NIF  designs  incorporate  very  thick  layers  of  DT,  which 
act  as  both  fuel  and  ablator  and  are  surroimded  by  a  very  thin  layer  of  normal- 
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density  polymer  material.  (The  polymer  is  required  to  contain  the  DT  vapor 
pressure.  It  is  unimportant  for  the  target  physics.) 

A  number  of  one-dimensional,  direct-drive,  point-design  calculations  have  been 
carried  out.  Figure  10  displays  the  results  for  four  values  of  a.  (The  value  of  a  is  the 
mass-weighted  average  over  the  cold  main  fuel  layer  at  a  time  close  to  ignition.)  In 
Fig.  10  tangential  focusing  is  assumed:  (i.e.,  the  foci  of  the  individual  laser  beams  are 
fixed  such  that  the  envelope  of  each  beam  covers  an  entire  hemisphere  of  the  initial 
capsule  surface). 

On  the  basis  of  the  simplified  model  presented  above,  a  conservative  approach  to 
designing  an  ignition-and-gain  direct-drive  capsule  design  for  the  NIF  would  begin 
with  an  a  between  3  and  4.  The  inner  and  outer  radii  of  the  0.25g/ cm^  DT  material 
are  1350  /an  and  1695  fim.  For  this  "all-DT"  design  there  is  an  outer  polymer  (CH) 
layer  with  thickness  of  4  /im.  The  initial  section  of  the  temporal  pulse  shape — ^the 
foot — ^has  been  designed  to  both  place  the  capsule  on  the  desired  implosion 
isentrope  and  minimize  the  development  of  the  Rayleigh-Taylor  instability  during 
this  foot  portion  of  the  implosion.  The  power  during  the  foot  rises  very  rapidly  from 
zero  to  a  constant  power  of  10  TW  and  is  held  constant  for  5  ns.  The  main  drive 
portion  is  temporeJly  shaped  to  achieve  the  desired  implosion  velocity  while 
minimizing  further  preheating  of  the  main  fuel  layer.  The  drive  portion  has  a 
temporal  duration  of  approximately  5  ns  with  the  major  portion  of  the  energy 
delivered  in  a  2.55  ns  section  with  a  peak  power  of  475  TW.  The  one-dimeirsional 
gain  calculated  for  this  design  is  ~20 
-  twice  that  quoted  for  indirect-drive 
capsule  designs  at  the  same  incident 
energy. 

The  results  presented  above  are 
from  cine-dimensional  simulations. 

The  next  step  is  to  address  the 
growth  of  the  perturbations.  A  large 
number  of  two-dimensional 
numerical  simulations  have  been 
carried  out,  using  one  perturbation 
wavelength  at  a  time  for  a  selection 
of  modes  £  in  the  range  of  10  to  1000. 

The  time-dependent  irradiation 
nonuniformity  in  each  mode  was 
obtained  from  a  (simplified)  two- 
dimensional  model  for  2-D  SSD  that 
agrees  with  more-detailed  three- 
dimensional  modeling  of  the  time- 


incident  laser  energy  (MJ) 

Fig.  10  Pellet  gain  versus  laser  energy  for 
the  all-DT  target,  for  four  values  of  the 
isentrope  parameter  a  =  P/Pp0 
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dependent  smoothing.  Typically  50  to  80  individual  perturbation  wavelengths  were 
modeled  for  each  capsule  design. 

The  evolution  of  these  modes  was  then  combined  with  a  saturation  model  to 
determine  the  temporal  evolution  of  the  unstable  mix  thickness  during  the  capsule 
implosion.  The  predicted  rear-surface  mode  spectrum  was  then  used  in  full  two- 
dimensional  simulations  to  determine  at  what  point  the  capsule  fails  to  ignite  and 
bum  (assuming  the  shell  did  not  totally  break  up  during  the  acceleration  phase). 
Figure  11  displays  the  results  of  these  calculations  for  the  a  =  3  design.  The 
simulations  showed  that  if  the  amplitude  of  the  inner  surface  perturbations  during 
the  deceleration  phase  exceeded  1.5  film  rms,  the  capsule  would  fail  to  ignite.  Using 
these  results,  one  can  arrive  at  a  "trade-off"  with  respect  to  the  levels  of  initial 
perturbations  the  a  =  3  design  can  tolerate  and  still  ignite  and  bum.  Using  the  1.5 
pm  iimer  surface  requirement  in  conjunction  with  modal  spectmms  measured  for 
the  outside-surface  initial-surface  finish,  and  the  two-dimensional  imprinting 
simulations,  one  finds  that  tihe  combined  perturbations  associated  with  initial- 
surface  nonuniformities  and  imprinting  must  be  maintained  below  1000  A  for  this 
capsule  to  still  ignite  and  bum  (a  value  very  close  to  that  obtained  using  the 
simplified  model).  A  "trade-off"  can  be  jjerformed  if  the  sources  of  pertiubations  are 
independent  and  therefore  add  in  quadrature.  If  one  asstimes  100  A  and  a  750  A 
initial  polymer  surface  finishes  then  the  level  of  imprinting  that  will  still  ignite  is 
respectively  995  A  and  660  A.  Direct  two-dimensional  numerical  simulations  of  the 
imprinting  with  2-D  SSD  on  this  capsule 
design  find  a  lower  imprint  level  of 
~  400  A  to  500  A  rms.  Thus  this  capsule 
design  does  have  the  potential  to  ignite 
and  biun  on  the  NIF.  The  two- 
dimensional  simulations  have  also 
foimd  that  at  these  levels  of  imprinting 
the  capsule  will  provide  essentially  ID 
yield  performance. 

This  target  design  has  also  been 
evaluated  independently  with  the  NRL 
ID  and  2D  implosion  codes.  Some 
differences  have  been  found  in  the 
timing  and  behavior  of  shocks  that 
bounce  off  the  origin,  but  the  target 
pulse  shape  can  be  readjusted  to 
produce  ignition  with  similar  laser 
nonuniformities. 


Fig.  11  Yield  versus  rms  perturbation  at 
the  inner  siurface  of  the  DT  ice/gas 
interface,  at  beginning  of  pellet 
deceleration,  showing  that  this 

perturbation  should  be  less  than  1.5  )jm 
for  ignition. 
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Some  improvements  in  both  target  gain  and  stability  have  been  investigated  at 
URLLE.  One  design  under  consideration  replaces  the  portion  of  the  DT  capsule  that 
is  ablated  by  the  laser  with  a  DT-wetted  polymer  foam  matrix.  This  results  in  an 
increase  in  gain  to  approximately  50  (due  to  higher  laser-energy  absorption)  with 
little  or  no  degradation  in  stability.  Additional  design  modifications  aimed  at 
reducing  the  levels  of  imprinting  are  being  actively  pursued. 

B.  Target  design  with  shock-preheated  foam  ablator 

1.  Pellet  Design  Parameters 

A  different  approach  to  controlling  the  RT  instability  uses  an  ablator  consisting 
of  low-density,  empty  CH  foam.  For  a  given  pressure  (from  the  first  shock)  the 
density  is  now  lower  and  the  isentrope  is  higher  compared  to  a  solid  density  CH 
ablator.  Shocks  moving  from  the  ablator  into  the  DT  will  partially  reflect  off  the 
interface  due  to  the  lower  shock  speed  in  the  DT.  This  reflected  shock  then  heats  and 
compresses  the  ablator  a  second  time.  For  example,  a  50  mg/cm^  foam  would  be 
shocked  a  factor  of  12  times  to  0.6  g/ cm^  after  the  first  shock,  while  tire  0.25  g/cm^ 
DT  fuel  would  increase  to  ~  1  g/cm^. 

A  possible  problem  with  this  approach  is  that  the  fuel/ ablator  interface  is 
classically  RT  unstable.  However  if  this  density  ratio  is  small  enough  after  the  iiritial 
compression,  then  the  RT  growth  may  be  acceptable.  For  the  above  situation,  the 
Atwood  number  is  (p2-pi)/ (P2+Pl)  =  0.4/1. 6  =  0.25.  With  this  low  density  jump  the 
2D  simulations  show  that  the  interface 
RT  irrstability  is  not  dominant. 

Figure  8b  illustrates  this  target 
design.  There  is  a  thin  vapor  barrier 
between  the  two  materials  that  has  been 
ignored  for  simplicity  in  the 
simulations.  The  studies  have 
emphasized  a  larger  laser  energy  of 
~  3.5  MJ,  to  examine  the  potential  of 
higher  energy  gains.  In  one  example,  the 
inner  radius  of  the  0.25  g/cm^  DT  is 
1800  pm.  The  irmer  and  outer  radii  of 
the  CH  foam  are  2077  pm  and  2947  pm. 

The  foot  of  the  pulse  rises  quickly  from 
zero  to  3  TW  (producing  2.75x1012 
•  W/cm2),  shocking  the  DT  to  about  2 
MBar.  The  foot  of  the  pulse  then  rises 
slowly  for  20  ns,  reaching  half  power  at 


Fig.  12  Density  structure  at  six  different 
times  during  the  implosion  of  the  target 
with  a  CH  foam  ablator. 
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26.4  ns  and  full  power  of  440  TW  at  28.1  ns,  dien  remaining  constant  for  7  ns. 

The  mass-averaged  shell  velocity  reaches  ~  3.3x10^  cm/sec  before  stagnation  and 
compression  to  460  g/cm^.  The  ID  gain  is  approximately  80.  Higher  gains  should  be 
possible  with  "zooming"  of  the  laser  beams.  (With  ISI,  one  can  easily  reduce  the 
laser  focal  spot  size  during  the  implosion  to  better  match  the  size  of  the  imploding 
pellet^.) 

2.  2D  Stability  Analysis 

This  target  design  was  analyzed  with  NRL's  FAST2D  Eulerian  sliding-zone 
implosion  code^l.  This  code  can  simultaneously  evaluate  several  perturbation 
modes  and  follow  them  into  the  ftonlinear  regime.  Figure  12  shows  one  of  the  2D 
calculations  of  laser  imprinting  and  RT  growth  using  a  2  THz  KrF  laser  with  ISI.  The 
bubbles  do  not  break  through  and  eventually  are  squeezed  out  due  to  the  spherical 
convergence.  The  inner  surface  of  tiie  shell  maintains  its  integrity  in  spite  of  the  ISI 
laser  imprinting.  Separate  calculations  show  that  the  target  cilso  maintains  its 
integrity  if  the  DT  inner  surface  roughness  is  less  than  0.5  jim  rms.  The  most 
sensitive  part  of  this  design  is  apparently  the  requirement  on  the  outer  surface  finish 
of  the  foam.  The  target  imploded  successfully  in  the  calculations  when  the  outer 
surface  perturbation  was  500  A  rms,  but  failed  when  the  surface  perturbation  was 
5000  A  rms.  The  best  foams  consist  of  short,  thin  fibers  that  are  separated  by  about 
1000  A,  and  will  introduce  very  high-mode  perturbations.  The  inherent 
perturbations  near  mode  number  100  are  unknown.  Target  fabrication  may  thus  be 
the  limiting  constraint  on  the  success  of  this  type  of  target. 

Figure  13  provides  a  useful  plot  of  the  minimum  and  angle-averaged  shell 
masses  per  steradian,  where  the  mass 
modulations  arise  from  laser  imprinting 
(time  moves  from  right  to  left  as  the 
shell  implodes  and  the  radius 
decreases).  The  shell  mass  in  the  bubble 
region  is  only  20%  less  than  the  angle- 
averaged  shell  mass.  Later  in  time  the 
bubbles  close  up  and  die  shell  mass  in 
the  bubble  region  approaches  the  angle- 
averaged  mass. 

This  implosion  was  recalculated  with 
a  bandwidth  of  0.3  THz,  with  all  other 
parameters  the  same.  The  inner  surface 
still  maintained  its  integrity,  although 
the  shell  mass  in  the  bubble  region 
dropped  to  40%  of  the  average  shell 


Fig.  13  For  the  target  with  a  CH  foam 
ablator,  average  (solid  line)  and 
minimum  (dashed  line)  mass  per 
steradian  during  the  implosion.  Time 
proceeds  from  right  to  left  as  the  pellet 
radius  decreases. 
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mass.  The  target  still  successfully  implodes,  although  one  would  expect  degraded 
performance. 

C.  Target  design  with  radiation-preheated  DT/CH-foam  ablator 

The  third  target  concept  uses  a  thin  layer  of  high-Z  material  such  as  tungsten 
surrounding  the  ablator.  Soft  x-rays  generated  during  the  foot  of  the  laser  pulse 
preheat  the  ablator.  The  ablator  can  not  be  pure  DT,  because  the  opacity  of  DT  is  too 
low  to  efficiently  stop  x-rays.  The  ablator  can  not  be  solid  CH,  because  the  opacity  is 
then  too  high,  and  the  soft  x-rays  would  only  penetrate  into  the  outer  part  of  the 
ablator.  We  therefore  propose  to  use  a  mixture  of  low-density  CH  foam  wicked  with 
frozen  DT.  The  ratio  of  CH  foam  mass  to  DT  mass  can  be  adjusted  to  optimize  the 
opacity  of  the  ablator.  For  sub  kilovolt  x-rays  the  minimum  opacity  of  carbon  occurs 
below  its  K-edge;  284  eV  for  cold  carbon,  (rising  to  ~  300  eV  at  higher  electron 
temperatures).  When  tungsten  is  heated  to  50-70  eV  by  the  foot  of  the  laser  pulse,  it 
is  an  excellent  emitter  of  radiation  near  250  eV,  with  ~  30  eV  width.  The  tungsten 
radiation  lasts  only  a  few  nanoseconds,  after  which  it  is  blown  off. 

The  detailed  design  of  this  target  hzis  had  to  await  the  development  of  improved 
non-LTE  opacity  models  that  can  evaluate  the  timgsten  radiation.  There  have  been 
several  recent  improvements  in  non-LTE  opacity  modeling^,  along  with  practical 
techniques  for  rapidly  calculating 
opacities  in  a  radiation-hydrodynamics 
code,  and  calculations  of  this  target 
have  now  begun.  The  optimum 
thickness  of  the  tungsten  is  in  the 
range  of  0.05-0.1  mg/cm^.  Its  radiation 
lasts  approximately  5  ns.  Several 
ablator  mixtures  have  been  considered 
of  the  form  CH(DT)n-  The  optimized 
opacity  is  in  the  vicinity  of  n  =  16-24,  so 
that  most  but  not  all  of  the  ablator  is 
heated  directly  by  the  soft  x-rays  from 
the  tungsten. 

Figure  14  plots  the  density  p  and 
isentrope  a  versus  radius  at  the  time 
when  the  series  of  shocks  first  break 
through  the  pellet  shell.  The  ablator 
material,  CH(DT)i6,  has  the  high 
isentrope  for  good  RT  stabilization, 
while  the  DT  fuel  is  on  a  very  low 
isentrope.  There  is  a  small  density 
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Fig.  14  Density  p  and  isentrope  a  profiles 
of  the  radiation  preheated  pellet  at  the 
time  when  the  shock  reaches  the  inner 
surface  of  the  shell,  demonstrating  that 
the  ablator  can  be  put  on  a  high  isentrope, 
while  the  DT  fuel  is  on  a  low  isentrope. 
The  laser  pulse  was  shaped  to  produce 
high  gain  target  performance. 
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jump  between  the  ablator  and  the  fuel,  due  to  the  radiation  preheating  of  the  ablator. 
FAST2D  calculations  indicate  that  the  laser  imprinting  is  a  factor  of  two  less  than  the 
imprinting  on  the  same  target  without  a  tungsten  coating.  Evaluation  and 
optimization  of  this  target  design  concept  is  ongoing  at  NRL  to  determine  if  a 
spatially  varying  isentrope  a  can  indeed  provide  sufficient  RT  stabilization  while 
simultaneously  maintaining  the  DT  fuel  on  a  very  low-isentrope. 

D.  Target  with  overcoats  of  foam  and  high-Z  shell. 

There  has  been  a  proposal  to  add  a  low-density  "foam  buffer"  and  a  high-Z  outer 
coat  to  the  outside  of  the  ablator  of  an  unspecified  pellet.  Prior  to  the  usual  foot  of 
the  laser  pulse,  there  would  be  a  separate  laser  pulse  that  heats  the  high-Z  material, 
producing  x-rays  that  propagate  supersonically  through  the  low  density  foam  and 
generate  a  long  scale  length  corona  on  die  ablator.  The  goal  here  is  to  reduce  the 
laser  imprinting.  The  use  of  supersonic  radiation  for  this  purpose  was  first  proposed 
by  M.  Desselberger  et  al.  in  Ref.  25.  More  general  suggestions  to  reduce  the  imprint 
with  x-ray  preheat  were  discussed  earlier  in  Refs.  22  and  67.  Numerical  simulations 
using  FAST2D  and  ORCHID  predict  that  an  x-ray  prepulse  reduces  but  does  not 
totally  eliminate  the  laser  imprinting.  Proponents  of  the  foam  overcoat  are  now 
attempting  to  incorporate  this  concept  into  a  self-consistent  ignition  target  design. 
(This  design  differs  from  the  one  in  Section  C  tiiat  uses  a  high-Z  overcoat  to  preheat 
a  DT-filled  CH  ablator  to  reduce  the  RT  instability.) 

IV.  LASER-PLASMA  INSTABILITIES 

Laser-plasma  interactions  can:  (1)  produce  suprathermal  electrons  that  preheat 
the  target  fuel  and  prevent  high  compression;  (2)  scatter  the  laser  light  away  from 
the  pellet,  and  reduce  the  coupling  efficiency;  (3)  filament  the  laser  light,  thereby 
increasing  the  peak  laser  intensity  and  enhancing  the  other  instabilities;  (4)  create 
denser  structures  in  the  xmderdense  plasma  that  decouple  the  laser  energy  from  the 
target. 

We  have  completed  an  initial  evaluation  of  the  three  major  decay  instabilities  for 
the  all-DT  and  all-CH  foam  ablator  designs  described  in  the  previous  Section.  We 
utilized  the  predicted  density,  temperature,  and  velocity  profiles  for  these  pellets  at 
various  times  during  the  implosion,  and  evaluated  how  far  the  designs  were  above 
or  below  the  effective  thresholds  using  models  that  assumed  a  coherent  laser  beam. 

Currently  three  important  elements  are  missing  from  the  theoretical  models. 

First,  the  total  laser  intensity  includes  all  of  the  overlapped  laser  beams  that  are 
incident  from  different  directions.  For  some  instabilities  we  would  expect  these 
different  laser  beams  to  behave  independently,  which  would  then  lower  the  effective 
laser  intensity  for  driving  laser-plasma  instabilities.  Second,  the  theories  assume  that 
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the  laser  has  a  constant  intensity.  However  these  beams  have  both  temporal  and 
spatial  fluctuations.  For  ISI  beams  the  speckle  satisfies  Gaussian  statistics^  (this  is 
also  nearly  true  for  broad-band  SSD,  except  for  the  lower  spatial  frequencies).  Under 
these  conditions,  the  probability  of  observing  a  given  intensity^^  is  exp(-Z/7o),  where 
lo  is  the  average  intensity.  For  incoherent  light,  4%  of  the  laser  energy  is  contained  in 
intensities  that  are  more  than  five  times  the  average  intensity.  If  this  4%  of  the  laser 
energy  is  converted  efficiently  into  fast  electrons  by,  for  instance,  the  two-plasmon 
decay  mode,  then  the  fuel  preheat  could  be  excessive.  (These  two  trends,  different 
incident  emgles  and  temporal  statistics,  will  tend  to  offset  each  other.)  Third,  decay 
instabilities  can  be  driven  out  of  resonance  when  the  laser's  frequency  moves 
outside  of  the  instability  resonance  width;  this  resonance  width  is  approximately  the 
growth  rate  y  For  some  of  the  decay  instabilities,  the  growth  rate  should  be 
reduced^o  from  y to  f-IA(o,  when  do  >  y  (Except  for  a  homogeneous  medium,  in 
which  the  reduction  of  growth  rate  is  compensated  by  an  increase  in  the  convective 
interaction  length^.)  These  three  physics  effects  are  in  need  of  further  theoretical 
and  experimental  research.  The  instability  analyses  in  the  following  subsections 
should  be  viewed  accordingly. 

A.  Filamentation 

A  laser  beam  can  filament  in  a  plasma,  with  an  enhanced  index  of  refraction 
induced  by  ponderomotive  forces  or  thermal  heating.  This  filamentation  instability 
can  be  seeded  at  a  finite  amplitude  by  speckle  in  the  laser  beams.  The  growth  of 
these  speckles  can  then  raise  the  peak  laser  intensities  and  enhance  the  various  laser- 
plasma  instabilities.  However  with  beam  smoothing  these  laser  speckles  are 
constantly  changing,  typically  on  a  time  scale  of  less  than  one  picosecond,  before 
they  can  undergo  significant  self-focusing.  Because  of  this  rapid  time  variation, 
along  with  a  plasma  corona  size  of  only  a  millimeter  (for  direct  drive  targets),  only 
modest  levels  of  time-dependence  are  sufficient  to  nearly  quench  the  filamentation 
instability.  Filamentation  should  not  be  important  for  direct  drive  targets,  using 
either  SSD  or  ISI.  There  are  theoreticaF2  and  numerous  experimental  studies^^,  74, 75, 
76, 77, 78, 79, 80  jjjat  support  this  conclusion. 

6.  Stimulated  Raman  Backscatter 

Stimulated  Raman  Scatter  (SRS)  is  the  decay  of  an  electromagnetic  wave  into  a 
plasma  wave  zmd  a  backscattered  electromagnetic  wave.  Assuming  a  linear  density 
profile,  Raman  backscattering  is  a  convective  instability  except  near  the  quarter 
critical  density  where  it  merges  with  the  two-plasmon  decay  mode  to  form  a  high 
frequency,  mixed  polarization,  absolute  instability®^.  For  densities  below  Uc/ 4,  the 
SRS  backscattering  threshold  for  significant  convective  gain®2/83  js  given  by 
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where  I  is  in  units  of  W/cm^,  L„  is  the  density  scale  length  in  units  of  100  p.m, 
and  Ao  is  the  laser  wavelength  in  units  of  jim. 

Using  this  formula,  the  two  direct-drive  targets  we  have  considered  are  below 
threshold  for  convective  SRS  growth  through  the  entire  under  dense  plasma,  at  all 
times  during  the  implosion,  because  of  the  relatively  large  density  gradients  and  the 
relatively  modest  peak  laser  power  with  direct-drive  targets.  This  does  not  exclude 
the  possibility  of  SRS  occurring  during  the  time-dependent  statistical  variations  of 
the  speckle  associated  with  beam  smoothing,  or  from  local  maxima  or  minima  (Ref. 
60)  induced  by  other  instabilities  such  as  the  RT  instability  or  the  two-plasmon 
decay  mode®^.  (Recent  impublished  OMEGA  pellet  experiments  without  beam 
smoothing  have  shown  negligible  SRS  when  several  beams  are  incident  from 
different  directions.) 

C.  Stimulated  Brillouin  Scatter 

Stimulated  Brillouin  Scatter  (SBS)  is  the  decay  of  an  electromagnetic  wave  into  an 
ion  acoustic  wave  and  a  backscattered  electromagnetic  wave.  In  practical  units  the 
SBS  threshold  for  significant  convective  growth  is®^/  83 


jsbs>i_7 


(n/njL^A^ 


(8) 


where  Lu  is  the  velocity  scale  length  in  imits  of  100  pm,  and  Tg  is  the  electron 
temperature  in  keV. 

This  formula  predicts  that  the  two  targets  considered  would  be  above  tihreshold, 
if  the  laser  intensity  is  the  stun  of  all  of  the  beams.  The  maximum  level  of  instability 
would  occur  as  the  laser  pulse  first  rises  from  the  foot  to  the  main  pulse;  this 
condition  has  a  low  temperature  and  a  high  laser  intensity.  For  the  all-DT  target 
with  1/3  pm  light,  I  / Ijgj  reached  a  value  of  6  at  a  density  of  0.2  Ug.  For  the  CH- 
foam  ablator  with  1  /4  pm  light,  I  /  reached  a  values  of  ~1  at  the  transition  to  the 
main  pulse. 

The  above  analysis  used  the  velocity  gradients  predicted  for  direct  drive  targets. 
Although  the  gradients  are  gentle  enough  so  that  the  instability  threshold  can  be 
exceeded,  the  gradients  are  also  large  enough  to  dominate  the  physics.  There  is  the 
possibility  that  the  temporal  variations  of  the  laser  speckle  pattern  may  drive 
acoustic  waves  that  stimulate  the  SBS  instability. 

Each  spot  on  the  surface  of  die  pellet  is  illiuninated  by  many  laser  beams,  which 
raises  the  question  of  power  transfer  between  beams.  This  effect  may  be  important 
for  other  types  of  laser  fusion  targets®^  that  require  beam  phasing  and  beam  overlap 
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at  very  high  laser  intensities.  However  for  direct-drive  this  power  transfer  should 
not  be  important  because  of  the  excellent  power  balance  between  beams,  and 
because  of  the  temporal  and  spatial  incoherence  between  the  beams. 

The  above  analysis  does  not  include  any  plasma  fluctuations  that  would  arise 
from  the  temporal  and  spatial  incoherence  in  the  laser  beams,  or  long  wavelength 
fluctuations  due  to  other  instabilities  such  as  RT  and  the  two  plasmon  decay,  which 
may  stabilize  SBS  (whereas  they  can  destabilize  SRS  as  noted  above). 

One  can  not  yet  conclude  that  there  will  or  will  not  be  significant  levels  of  SBS 
instability  in  these  pellets.  This  will  be  a  subject  for  further  experimental  and 
theoretical  research. 

D.  Two-Plasmon  Decay 

The  two-plasmon  instability  is  the  decay  from  an  electromagnetic  wave  into  two 
electron  plasma  waves,  and  occurs  near  quarter  critical  density^^.  In  practical  units, 
the  threshold  condition  is 


>0.54^^  (9) 

With  the  above  formula,  this  instability  is  potentially  very  dangerous  for  both 
pellet  designs.  However  the  peak  growth  rate  of  this  instability  for  short  wavelength 
lasers  has  a  long  transverse  wavelength,  and  it  may  be  quenched  by  th.e  spatial 
incoherence  in  the  laser  beams.  Given  that  the  threshold  is  so  low,  many  modes  will 
be  driven  unstable  and  the  most  dangerous  modes  may  be  those  that  are  not  most 
resonant,  but  which  suffer  large  Landau  damping  and  extend  to  densities  well 
below  nc/4.  Such  modes  will  no  longer  preserve  the  geometry  of  the  most  unstable, 
highly  resonant  modes,  but  can  slowly  grow  and  convectively  satiurate  at  high  levels 
long  after  the  more  resonant  modes  have  saturated  nonlinearly.  However  these  non 
resonant  modes  with  low  growth  rate  may  be  more  easUy  stabilized  by  the  broad 
bandwidth  of  lasers  with  ISI  or  SSD. 

The  above  analysis  is  supported  by  recent  high-intensity  OMEGA  pellet 
experiments  that  used  imsmoothed  laser  beams.  There  was  no  evidence  for  the 
resonant  two-plasmon  mode  when  multiple  beams  were  overlapped,  but  there  was 
evidence  of  the  non  resonant  modes. 

This  instability  has  the  potential  of  generating  harmful  energetic  electrons 
tiirough  Landau  damping  and  wave  breaking  of  plasma  waves.  Neither  the  details 
of  this  instability  nor  the  electron  generation  process  are  well  understood.  Further 
studies  are  necessary,  to  include  the  effects  of  multiple  laser  beam  overlap  with 
broad  bandwidth  light,  and  the  probabilistic  variation  in  laser  intensity. 

E.  Radiating  Plasma  Structures 
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If  a  target  contains  material  other  than  DT,  then  radiation  can  have  two  effects. 
On  the  positive  side,  it  can  reduce  the  density  gradients  and  reduce  the  RT  growth 
rate.  On  the  negative  side,  computational  modeling  with  highly-resolved  opacity 
data  sometimes  predicts  the  occurrences^  of  Radiating  Plasma  Structures  (RPS). 
These  RPS  are  related  to  the  radiative  condensation  instabilitySS.  They  occur  when 
the  laser  intensity  is  a  function  of  time,  and  the  function  L  =  ApK^cfT*  goes  through  a 
maximum  as  a  function  of  T,  where  a  is  the  Stefan-Boltzmann  constant,  p  is  the 
density,  and  Kp  is  the  mean  Planck  opacity.  Then  plasmas  at  a  lower  temperature  can 
radiate  more  than  adjacent  plasmas  at  a  higher  temperature.  The  lower  temperature 
material  forms  a  local  higher  density  region  that  can  then  interfere  with  the  laser 
propagation.  For  carbon,  the  main  effect  occurs  when  the  electron  temperature  is  in 
the  vicinity  of  65  eV,  and  the  K-shell  of  carbon  is  excited.  This  RPS  can  complicate 
the  use  of  some  otherwise  desirable  laser  pulse  shapes.  The  structures  are  sometimes 
seen  in  radiation-hydrod5mamics  codes,  but  have  not  yet  been  directly  observed  in 
experiments. 

F.  Fokker-Planck  Modifications  to  Thermal  Transport 

The  computer  models  that  are  used  to  design  targets  assume  that  thermal 
transport  satisfies  a  Spitzer-Harm  diffusion  equation.  However  it  has  been  noted®^/ 
90, 91,  92, 93,  94  that  even  if  the  electron  distribution  function  is  a  Maxwellian  in  die 
absorption  region,  the  electrons  in  the  tail  can  have  a  mean  free  path  that  reaches 
farther  than  the  thermal  diffusion  front,  producing  preheat  perhaps  all  the  way  into 
the  DT  fuel.  There  has  been  a  recent  attempt  to  model  these  kinetic  effects  on 
ignition  target  design^®,  but  improvements  will  be  needed  in  the  modeling.  We 
therefore  can  not  yet  quantitatively  evaluate  the  potentially  deleterious  effects  of  this 
fast  electron  preheat. 

V.  LASERS  FOR  A  REACTOR 

For  a  fusion  reactor  using  direct-drive  pellets,  die  two  lasers  currently  imder 
consideration  are  a  KrF  laser  and  a  DPSSL.  These  lasers  will  have  to  satisfy  severe 
requirements  in  efficiency,  repetition  rate,  cost,  durability,  and  laser  beam  quality. 
Neither  system  can  now  meet  all  of  the  requirements,  but  the  advocates  of  the  two 
lasers  have  proposed  solutions  that  may  make  them  viable.  Because  of  the  modular 
design  of  fusion  lasers,  consisting  of  many  parallel  beam  lines,  it  should  be  possible 
to  develop  and  evaluate  both  laser  concepts  with  modest  levels  of  funding. 

As  an  example,  a  KrF  driver  could  utilize  100  large  amplifiers,  each  producing  34 
kj  of  laser  energy.  Each  of  these  amplifiers  in  turn  would  be  pumped  by  four 
identical  electron  beam  systems,  each  producing  an  "equivalent"  of  8.5  kJ.  This  core 
module  would  be  only  slightly  larger  than  the  5-6  kJ  achieved  with  the  Nike  laser  at 
NRL.  The  existing  pulsed  power  technology  would  be  replaced  by  efficient  and 
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reliable  non-linear  magnetic  switches  similar  to  the  Repetitive  High  Energy  Pulsed 
Power  device^^  (RHEPP),  or  the  recently  developed  Semiconductor  Opening 
Switch^^  (SOS).  The  optical  design  of  a  KrF  fusion  laser  could  be  similar  to  Nike,  and 
it  would  use  the  same  type  of  gas  mixture  of  Ar/Kr/F2- 

For  the  DPSSL,  the  Xe  flash  lamps  used  in  lasers  such  as  OMEGA  and  the  NIF 
would  be  replaced  by  efficient  and  reliable  semiconductor  diode  arrays^^.  The 
Ndrglass  lasing  medium  would  be  replaced  by  a  rare-earth  doped  crystal  with  a 
longer  lived  upper  state  and  an  absorption  band  chosen  to  effectively  use  the  diode 
arrays.  Turbulent  gas  cooling  has  been  proposed  to  cool  the  laser  crystals.  The  laser 
architecture  would  use  the  same  type  of  multi-pass  extraction  as  the  NIF. 

A.  Efficiency 

The  reactor  study  called  Sombrero^^  estimated  the  total  efficiency  of  a  KrF  laser 
to  be  7.5%,  plus  2%  heat  recovery  from  the  laser  gas.  This  requires  a  75%  efficient 
pulsed  power  system,  compared  to  the  existing  60%  of  RHEPP  and  50%  of  SOS. 

Reference  98  presented  a  design  for  a  DPSSL  with  a  total  laser  efficiency  of  8.6%. 
This  study  assumed  a  diode  efficiency  of  55%,  improved  from  the  level  of  current 
diodesi^  of  45-50%.  There  is  promising  recent  research^oi  that  has  ridsed  the  diode 
efficiencies  to  as  high  as  66%. 

B.  Repetition  Rate 

The  diode  pump  should  easily  run  at  a  high  repetition  rate.  In  fact,  a  higher 
repetition  rate  increases  the  duty  cycle  of  the  diodes,  and  thus  reduces  their  capital 
cost  in  comparison  to  the  fusion  power  output.  However  the  beam  quality  of  a 
DPSSL  xmder  appropriate  loading  must  still  be  proven.  This  will  be  evaluated  as 
part  of  the  new  Mercury  project  at  LLNL  (Lawrence  Livermore  National  Laboratory) 
that  utilizes  face  pumping  and  face  cooling  with  turbulent  helium-gas  flow^o^. 

For  KrF,  the  RHEPP  runs  at  120  Hz,  far  above  tiie  5-7  Hz  needed  for  a  reactor.  It 
uses  all  solid  state  components  emd  should  also  have  an  extremely  long  lifetime 
(>10^  shots).  The  SOS  has  a  repetition  rate  of  150  FIz.  The  lifetime  is  reported  to  be 
very  long.  The  EMRLD  electron-beam  pumped  excimer  laser  has  already 
demonstratedi03  amplification  of  a  laser  beam  at  only  1.3  times  diffraction  limited  at 
100  Hz. 

C.  Cost 

The  current  cost  of  KrF  power  supplies  would  be  dominated  by  the  high  price  of 
Metglass  that  is  used  in  the  saturable  inductors.  At  the  current  cost  of  $100/kg,  the 
cost  of  a  KrF  laser  would  be  too  high  by  a  factor  of  about  40.  However  Metglass  is 
slowly  being  adopted  by  the  power  industry,  and  costs  are  predicted  to  drop 
significantly.  (The  same  assumption  is  used  in  the  design  of  heavy  ion  accelerators 
for  a  fusion  reactor.) 
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The  current  cost  of  laser  diodes  is  in  excess  of  10  $/W,  too  high  by  a  factor  of  100. 
The  cost  is  currently  high  because  production  rates  are  low.  Costs  should 
dramatically  reduce  as  demand  increases  and  processing  steps  become  more 
automated.  Newer  diode  arrays  are  being  developed  that  are  more  compatible  with 
automation.  Scaling  laws  predict  costs  of  0.10  $/ peak-Watt  with  production  levels 
appropriate  to  one  reactor  per  year.  If  costs  drop  well  below  0.10  $/peak-Watt,  then 
one  could  use  more  diodes  per  laser  disk,  and  pump  the  laser  for  a  shorter  period  of 
time.  There  would  then  be  less  spontaneous  decay  from  the  excited  state,  and  a 
resulting  higher  laser  efficiency.  Efficiencies  as  high  as  15%  are  believed  possible  if 
the  cost  of  the  diodes  becomes  negligible  compared  to  the  other  reactor  components. 

D.  Durability 

A  very  attractive  feature  of  the  DPSSL  driver  is  its  durability.  Major  components 
are  expected  to  last  >  10^  shots,  except  for  routine  repolishing  and  replacement  of 
optical  components. 

KrF  lasers  aire  now  durable  enough  for  laser  fusion  resecirch,  but  not  for  a  power 
plant.  The  two  major  durability  problems  with  the  Nike  laser  are  the  lifetime  of  the 
thin  foils  that  separate  the  laser  gas  cell  from  the  electron  beam  source,  and  the 
lifetime  of  the  anti-reflection  coating  on  the  windows  of  tihe  laser  cell.  However 
longer  foil  lifetimes  were  achieved  on  the  EMRLD  laser;  their  foils  did  not  degrade^ 
even  after  10^  shots.  The  anti  reflection-coated  laser  windows  operate  in  a  rather 
hostile  environment  of  fluorine,  intense  UV  laser  light,  x-rays,  and  energetic 
electrons.  The  coatings  on  the  fused  silica  windows  of  the  large  Nike  amplifier 
degrade  after  a  few  hundred  laser  shots.  We  do  not  yet  know  the  reason  for  this 
degradation.  Solutions  have  been  proposed,  such  as  the  use  of  fluorine  resistant 
CaF2  windows,  and  it  is  the  subject  of  low-level  ongoing  research. 

E.  Laser  Beam  Uniformity 

Successful  pellet  implosions  require  optical  beam  smoothing  in  both  the  high 
transverse  modes  and  the  low  transverse  modes.  A  KrF  laser  with  ISI  provides  the 
maximum  level  of  beam  smoothing  in  botih.  regimes,  because  it  can  have  a  large 
bandwidth,  up  to  3  THz,  and  because  the  ''beamlets"  in  ISI  are  all  statistically 
independent.  Nike  has  already  demonstrated  the  high-mode  beam  uniformity  that 
would  be  necessary  for  a  fusion  reactor,  with  1%  rms  fluctuations  in  one  beam  in  a  4 
ns  pulse  at  1  THz.  This  extrapolates  to  ~  0.15%  rms  with  39  overlapping  beams. 
There  has  also  been  occasional  operation  at  3  THz. 

The  proposed  solid  state  lasing  medium,  Yb^^-  doped  Sr5(P04)3F,  (or  Yb:S-FAP), 
has  a  narrow  line  width,  producing  only  -0.3  THz  at  the  third  harmonic.  This  crystal 
was  chosen  to  maximize  the  diode's  pump  duration  while  still  producing  adequate 
gain'®'^.  There  are  other  possible  hosts^®^  for  Yb,  including  Ca  substituted  Yb;S-FAP 
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that  would  have  three  times  the  bandwidth.  However  the  Einstein  relationship 
shows  that  a  higher  bandwidth  implies  a  higher  saturation  fluence.  Changing  to  this 
higher  bandwidth  host  would  increase  the  saturation  fluence  from  3.2  J/cm^  to  6.0 
J/ cm^.  Since  efficient  extraction  of  the  inversion  requires  a  laser  fluence  that  is  about 
three  times  the  saturation  fluence,  this  would  require  higher  damage  threshold 
optical  components. 

The  beam  uniformity  of  2-D  SSD  has  so  far  been  limited  by  the  bandwidth  of  the 
frequency  converter,  and  the  use  of  a  single  color  cycle,  which  does  not  provide 
temporal  smoothing  of  the  lowest  mode  numbers.  The  frequency  conversion 
bandwidth  may  be  increased  by  using  advanced  two-tripler  schemes^^,  and  the  low¬ 
mode  distortion  may  be  reduced  by  increasing  the  munber  of  color  cycles  across  the 
beam^^.  The  potential  of  advanced  two-tripler  schemes  and  multiple  color  cycling  is 
under  investigation. 

VI.  PROSPECTS  FOR  FUSION 

Solutions  have  been  developed  to  the  problems  in  direct-drive  of  laser 
imiformity,  hydrodynamic  instability,  and  laser-plasma  instabilities.  There  are  now 
tiiree  specific  target  designs  that  may  meet  all  of  the  physics  constraints  and  achieve 
ignition.  The  target  design  codes  have  extensive  physics  packages,  and  they  have 
been  successfully  compared  with  numerous  laser-target  experiments  using  a  variety 
of  sophisticated  diagnostics.  The  remaining  disagreements  may  be  attributable  to  3D 
effects  that  are  not  in  the  2D  codes,  and  that  probably  require  that  the  2D  target 
desigiis  contain  a  margin  of  safety  for  larger-than  predicted  RT  saturation  effects. 
Two  of  the  target  designs  keep  the  DT  fuel  on  a  low  enough  isentrope  so  that  gains 
above  100  may  be  possible.  All  three  target  designs  may  be  applicable  to  the  NIF 
laser  that  is  now  under  construction. 

Although  the  prospects  for  ignition  with  direct  drive  are  now  very  promising, 
there  are  also  some  remaining  imcertainties.  The  Rayleigh-Taylor  instability  needs  to 
be  further  evaluated  in  plcinar  geometry  with  the  specific  ablator  materials  that 
would  be  used  in  ignition  targets.  The  Rayleigh-Taylor  instability  should  also  be 
evaluated  in  spherical  geometry,  to  evaluate  higher  laser  intensities  and  lower-mode 
perturbations.  Laser-plasma  instabilities  need  re-evaluation  with  statistically  smooth 
laser  beams  to  examine  the  roles  of  laser  beam  overlap  and  laser  bandwidth.  A 
larger  data  base  is  needed  to  better  evaluate  the  new  non-LTE  models  of  opacity  and 
emissivity.  The  design  of  the  three  targets  is  still  incomplete,  with  research  especially 
needed  into  Fokker-Planck  effects,  3D  RT  saturation,  final  deceleration  and  burn, 
and  non-LTE  modeling  of  emissivity.  Development  has  begun  on  a  rep-rate  DPSSL, 
and  a  parallel  program  should  be  started  to  develop  a  rep-rate  KrF  laser. 

The  different  capabilities  of  Nike  and  OMEGA  complement  each  other.  The  Nike 
KrF  laser  has  the  most  uniform  laser  beam,  and  was  optimized  for  a  long  pulse 


duration  and  modest  laser  intensity  to  study  laser  imprinting  and  the  Rayleigh- 
Taylor  instability  in  planar  geometry  scaled  from  fusion  targets.  The  OMEGA 
Nd:glass  laser  with  3a)  frequency  conversion  is  the  first  fusion  laser  with  a  sufficient 
number  of  laser  beams  to  properly  illuminate  and  symmetrically  implode  a  spherical 
pellet,  and  it  can  more  easily  reach  the  higher  laser  intensities  needed  for  the  study 
of  laser-plasma  instabilities. 
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